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In this work, we propose a new and simple model that supports Chern semimetals. These new 
gapless topological phases share several properties with the Chern insulators like a well-defined Chern 
number associated to each band, topologically protected edge states and topological phase transitions 
that occur when the bands touch each, with linear dispersion around the contact points. The tight- 
binding model, defined on the Lieb lattice with intra-unit-cell and suitable nearest-neighbor hopping 
terms between three different species of spinless fermions, supports a single Dirac-like point. The 
dispersion relation around this point is fully relativistic and the 3x3 matrices in the corresponding 
effective Hamiltonian satisfy the Duffin-Kemmer-Petiau algebra. We show the robustness of the 
topologically protected edge states by employing the entanglement spectrum. Moreover, we prove 
that the Chern number of the lowest band is robust with respect to weak disorder. For its simplicity, 
our model can be naturally implemented in real physical systems like cold atoms in optical lattices. 
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Introduction:- Topological phases represent one of 
most exciting and interesting field of condensed matter 
physics. Topological insulators and superconductors are 
well-known examples of free fermion systems defined by a 
gapped bulk that support robust gapless edge states [T]. 
All these systems described by non-interacting Hamilto¬ 
nians fit in the periodic table of fermionic topological 
phases and the topological phase transitions occur when 
the bulk gap closes. In particular, two-dimensional chi¬ 
ral topological systems like the Haldane model [2] and 
p-wave superconductors [3] are characterized by the topo¬ 
logical Chern number that fixes the number of topologi¬ 
cally protected edge modes. The former is an example of 
topological insulator in the class A, called also Chern in¬ 
sulator, where both time-reversal and particle-hole sym¬ 
metries are broken while the latter is a topological su¬ 
perconductor living in the class D, where particle-hole 
symmetry is instead preserved. 

In the last years, part of the research of new topological 
phases has focused on gapless bulk systems like three- 
dimensional Weyl [4] and Dirac semimetals [5] that rep¬ 
resent three-dimensional versions of graphene. In these 
systems, the bands touch each other in a discrete set 
of points that can be seen as point-like defects in mo¬ 
mentum space. These Weyl and Dirac points are topo¬ 
logically protected by well-defined Berry phases and the 
boundaries support suitable gapless modes. Clearly, 
when the time-reversal symmetry is broken in two di¬ 
mensions, the Chern number can characterize the bands 
of the semimetals only when those bands do not touch 
each other, i.e. when there are no Dirac or Weyl cones. 
Some extended Haldane models with these characteris¬ 
tics have been recently analyzed [SIE]. 

The goal of this work is to present a new and 
simple model that supports two-dimensional Chern 
semimetallic phases, which share several properties 



FIG. 1: (Up) This picture shows a square lattice with 

three species of fermions represented by triangles, circles and 
squares. The dashed circle encloses the unite cell. Intra-unit¬ 
cell hopping terms are represented by red lines, the nearest- 
neighbor ones with blue lines. The dashed green line identifies 
the hopping between non-adjacent fermions. (Low-left) Dis¬ 
persion energy for K — —J — 1 and M = 0. (Low-right) 
Projection of dispersion energy for with M = 0.5z. 


with the Chern insulators. Our tight-binding model is 
defined on a Lieb lattice with intra-unit-cell and suitable 
nearest-neighbor hopping terms between three species 
of free spinless fermions. The model supports only a 
Dirac-like cone due to the presence of a zero-energy 
flat band in the middle. In fact, it is possible to avoid 
the femion doubling in the lattice if a flat band is 
present. At this point, it is possible to deform the 
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bands by introducing a further hopping term between 
non-adjacent fermions in the unit cell. In this case, as we 
show in the next section, the lower band is characterized 
by a non-zero Chern number z/ = ±1 and the model 
supports robust edge states. We analyze the edge 
modes by the employing entanglement spectrum and we 
demonstrate the robustness of the Chern number with 
respect to the presence of weak disorder. Moreover, we 
show that the corresponding effective Hamiltonian heff 
is fully relativistic but different with respect to a Dirac 
Hamiltonian, because the 3x3 matrices in hgff satisfies 
the Dufhn-Kemmer-Petiau algebra [8] . This implies that 
our model does not fit in any already known periodic 
table of topological gapless phases based on K-theory 
and Clifford algebra [9]. Finally, due to the presence 
of only intra-unit-cell and suitable nearest-neighbor 
hopping terms, our model can be easily implemented 
in real physical systems like cold atoms in optical lattices. 

Lattice model and Chern number:- To begin with, 
we introduce our tight-binding model. We consider three 
species of spinless fermions on a Lieb lattice as shown in 
Fig. 1, described by the following Hamiltonian 

H = ^ ^ + b^^Ci)-\rK[a^^bi-^x + ^i^i+y) T 

i 

+h.c., (1) 

where i is the site index, x = (1,0), y = (0,1), and 
a, b and c are the three species of fermions represented 
in Fig. 1 in terms of triangles, circles and squares, re¬ 
spectively. Here, the tunneling coefficients J and K are 
taken real, while M = is complex. By imposing 

periodic boundary conditions we introduce the Fourier 
transformation where /c = 1,2,3 is 

the species index, to obtain H = Xlpwhere 
ipp = {ap^bp^CpY"^ p G BZ = [—7r,7r) x [—7r,7r) and the 
kernel h{p) is a 3 x 3 hermitian matrix. It is straightfor¬ 
ward to see that the time-reversal symmetry is broken, 
namely h{p) ^ h'^{—p), but is restored when M is real. 
In particular, when M = 0, the system supports a single 
Dirac-like cone and a zero-energy flat band as shown in 
Fig. 1. First of all, the presence of the flat band allows 
to avoid the fermion doubling problem as shown in m- 
In this case, our model behaves similarly to other three- 
band models define on Lieb and Kagome lattices m 
ng. However, a crucial difference emerges with respect 
to the latter, when a non-zero complex M is switched 
on. In this case, the complex hopping term with a small 
m < 1 opens a pseudo-gap between the lower and upper 
bands and deforms the middle bands such that the sys¬ 
tems remains gapless as shown in Fig. 1. We now show 
that in this semimetallic phase, our model behaves as a 
non-trivial topological semimetal, where each band has 
a well-defined Chern number z/. This is indeed possible 
because the time-reversal symmetry is broken and the 



FIG. 2: Energy bands for X = — J = 1 on cylinder of length 
Ly = 80 as a function of momentum px- (Up-left) M = 0. 
A single Dirac-like cone and a zero-flat band appear in the 
system. (Up-right) m — 0.5 and 0 — 7r/2. Upper and lower 
bands have non-zero Chern number and the crossing red lines 
represent the edge modes between them. (Low-left)?TZ =1.3 
and 0 — 7r/4. the upper gap opens and the Chern number is 
acquired by the middle band after a topological phase transi¬ 
tion when the bands touch. (Low-right) When even the lower 
gap opens, the system becomes a trivial insulator. 

bands do not touch each other at any point in the BZ 
like for example in graphene and Weyl semimetals. The 
Chern number related to the n-th band with normalized 
Bloch wave function \n{p)), is defined by 

— X : [ d pFxy^ (2) 

where the Berry connection (<a = x,y) and the cor¬ 
responding curvature tensor Fxy are given by = 
(n(p)|g|^|n(p)) and F^y = d^Ay - dyA^. Here, we are 
assuming that there is no degeneracy for the n-th state 
and we obtain a Chern number z/ = ±1 by integrating the 
Berry curvature on the lower band. For our numerical 
calculations, we use the discretised version of (§ derived 
in m- From a geometric point of view, the presence 
of the Chern number is connected to the presence of a 
non-zero flux passing through the triangles defined by 
the hoppings J and M inside the unit cell when M is 
complex. Fig. 3 reproduces the phase diagram of this 
new topological phase for K = 1 and 0 = nj2 and dis¬ 
plays how the chirality depends on the sign of m and 
J. Importantly, topological phase transitions occur when 
the bands touch each other, with linear dispersion around 
the contact points. Moreover, as we show in the next sec¬ 
tion, topologically edge states are associated to the non¬ 
zero Chern number. Thus, this topological semimetal 
shares several properties with the Haldane model but 
there are also several differences. Firstly, in the momen¬ 
tum space, there appears only a Dirac-like point. This 
property is due to the presence of a zero-energy fiat band 
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FIG. 3: (Up) Chern number phase diagram in function of 
J and m with K = 1 and 0 = tv j 2 when the lowest band is 
completely filled. (Low) Phase diagram in function of m and 
6>withK = -J=l. 

for M = 0 and is shared with some three-band realiza¬ 
tions of Chern insulators, even if in these last cases, the 
topological phases are obtained only by introducing spin- 
orbit interactions mHH]. Secondly, in the low-energy 
regime, the corresponding effective kernel Hamiltonian 
/left is not defined by any Dirac or Dirac-Weyl Hamil¬ 
tonian even if there is a single Dirac-like point. This 
is possible because /leff is a fully relativistic first-order 
Hamiltonian, given by 

/left = K [13^, l3°]p, + K[/3y,l3°]py + M 13°, (3) 

where we have fixed J = —K. Here, the 3x3 matrices 

/ 00 -l\ / 000 \ 

/3°= 00 0 , ^^= 00i , 

\100/ \0 -i 0 J 

/O -i 0\ 

= h 0 0 , (4) 

\ 0 0 0 / 

satisfy the following conditions 

where 77 ^^ is the the relativistic Minkowski metric such 
that diag 77 ^^ = (— 1 , 1 , 1 ) and there is no summation on 
repeated indexes. The above relations identify a gen¬ 
eralized Clifford algebra, called Duffin-Kemmer-Petiau 
algebra [ 8 ]. This algebra is associated to the Duffin- 
Kemmer-Petiau theory that describes relativistic spin-0 
and spin -1 particles by employing the same formalism 


used by Dirac for spin- 1/2 particles. In particular, the 
effective Hamiltonian in (§ formally describes spin -0 
quasiparticles in two dimensions and the corresponding 
spinor field satisfies the Klein-Gordon equations. This 
also implies that our model does not fit in any already 
known periodic table of topological semimetals because 
this classification is based on the K-theory associated to 
the standard Clifford algebra [9]. 

Edge states and entanglement spectrum:- As 
is well-known, one of the main properties of topologi¬ 
cal phases concerns the existence of robust edge states. 
In particular, in two-dimensional time-reversal broken 
phases, the number of chiral edge modes coincides with 
the value of Chern number. This is well established when 
the systems have a gapped bulk but still an open ques¬ 
tion for Chern semimetals. In order to show the edge 
state energy dispersion we impose open boundary con¬ 
ditions in the 7/-direction and periodic in the x-direction 
we perform a Fourier transformation that decomposes 
our two-dimensional Hamiltonian into decoupled one¬ 
dimensional Hamiltonians describing chains of length Ly 
parametrised by Px- In Fig. 2 it is clear that edge states 
appear and cross between bands with non-trivial Chern 
numbers of opposite sign. Here, we show that Chern 
semimetals have topologically protected edge modes by 
employing the entanglement spectrum m- Let us briefly 



FIG. 4: (left) The gap-crossing line for K = — J — l^m — 0.5 
and 0 = 7r/2 shows the presence of an edge mode in the 
entanglement spectrum on the cylinder, that represents the 
signature of the physical edge state, (right) The absence of 
the crossing line for K — —J = 1, tti ^ 1 and 0 — 7r/2 is 
due to the absence of edge states in the topologically trivial 
insulating phase. 

summarize the main properties of the entanglement spec¬ 
trum for free-fermion systems. First of all, we divide 
our model which is placed on a cylinder in two non¬ 
overlapping sub-regions through a cut that runs along 
X, the periodic dimension m- The reduced density ma¬ 
trix of the ground state p of the sub-system is related to 
the entanglement Hamiltonian i.e. p = and the 
corresponding eigenvalues contain the main informa¬ 
tion about topological phases. However, in the case of 
free fermion Hamiltonians, it is possible to prove that 
are in one-to-one correspondence with the eigenvalues Xi 
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of the correlation matrix m 

Cnm ~ (^'^| 5 (5) 

where \gs) is the ground state, the labels m^n are 
restricted to be within the sub-system, and in our case, 
hn represent the fermion operators a^, bn and c^. Thus, 
the entanglement spectrum in the model is the spectrum 
of the correlation matrix. Since our cut runs along the 
periodic direction x we plot the eigenvalues X{px) of the 
correlation matrix C{px) of each of the decoupled chains 
as showed in Fig. 4, where the ground state used to 
construct the correlation matrix is built by filling the 
lower band, and the entanglement cut is placed in the 
middle of each chain. In this figure, the virtual mode 
in the entanglement spectrum in the plot on the left 
represented by a gap-crossing line is the signature of 
the presence of a physical dispersing edge mode that 
will appear in the energy spectrum when a physical 
boundary is introduced in the place of the entanglement 
cut. This edge state is associated to the completely filled 
lower band which shows u = 1. Instead, on the plot on 
the right, clearly the gap-crossing line disappears in the 
topologically trivial insulating phase when the band has 
z/ = 0 . 

Disorder:- In this section, we show that the Chern 
semimetals are even robust phases in presence of disorder, 
which would at least weakly appear in the case of experi¬ 
mentally realising the model with cold atoms in an optical 
lattice. We impose periodic boundary conditions and de- 
vide the systems into four regions: three non-overlapping 
regions X, T, Z with common boundaries that form a 
tripple point which are all surrounded by the fourth re¬ 
gion. We then numerically calculate the Chern number 
in real space [201 ED 

u = 127ri EEE {CjkCkiCij - CjiCikCkj) (6) 

jex keY lez 

which gives the same value as equation in the trans- 
lationally invariant case. Numerically we are restricted 
to finite system sizes. However we use system size 
scaling in order to see that the real space Chern number 
converges to an integer as the system size is increased 
and thus choose an acceptable system size in order to 
perform our computation. Disorder is introduced in 
all of the hopping amplitudes T = X, J, m in the form 
T{1^ Jd) where Jd is the disorder amplitude. For each 
value of Jd we compute the Chern number for a number 
of disorder realisations and average u over them. The 
result as a function of 0 and Jd is shown in Fig. 5 
(up). For weak disorder, Jd < ^ the topological phase is 
robust showing non-trivial value of the Chern number. 
Since there is no energy gap in the semimetallic phase, 
we look at the density of states in order to identify 
the pseudo-gaps, where we have averaged a number 



FIG. 5: (Up) Real space Chern number u against 0 and 
Jd for K = —J — 1 and m — 0.5. Here a system size of 
20 X 20 unit cells is chosen and the particle filling is what 
corresponds to filling the lowest band in the clean case. Each 
disorder realisation was averaged 20 times, enough to discern 
the topological regions from the trivial ones. For weak dis¬ 
order Jd < 0.5 the topological phase is robust and for strong 
disorder the Chern number vanishes. (Low) Density of states 
p{E) for disorder amplitudes 0 < < 0.5 (inset) where we 

averaged 10^ disorder realisations on 26 x 26 unit cells for each 
Jd- The two minima correspond to the semimetalicl pseudo¬ 
gap. The density of states profile does not change around the 
pseudo-gaps. Here K = — J = 1, m = 0.5 and 0 = 7r/2. 


of disorder realisations for each disorder amplitude 
Fig. 5 (low). In the clean case the density of states 
shows minima at the energies where the pseudo-gaps are 
located. The minima remain for weak disorder and we 
only notice a spread in the extreme energies away from 
the pseudo-gaps. This behaviour agrees with the Chern 
number remaining \v\ « 1 for weak disorder showing the 
robustness of the Chern semimetallic phase. 

Conclusions:- Summarizing, we have shown that a 
new and simple tight-binding model on the Lieb lattice 
with only intra-unit-cell and nearest-neighbor hopping 
terms supports Chern semimetallic phases. These new 
topological semimetals are characterized by a non-zero 
Chern number in the bulk and topologically protected 
gapless edge states. We have proved the existence of 
the latter by studying the entanglement spectrum on 
the cylinder. Moreover, we have shown that the Chern 
semimetal is robust also in presence of weak disorder 
by calculating the Chern number in the real space. 
Finally, for its simplicity, our Lieb lattice model can be 
experimentally realized with cold atoms in an optical 
lattice [HEg. 
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